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Invariants for Darboux transformations 
of Arbitrary Order for 



D x D y + aD x + bD y + c 

Q-r 

i-Q 1 Ekaterina Shemyakova 

Abstract. We develop the method of regularized moving frames of Fels 
and Olver to obtain explicit general formulas for the basis invariants 
that generate all the joint differential invariants, under gauge transfor- 
mations, for the operators 
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■ C = D x D y + a(x, y)D x + b(x, y)D y + c(x, y) 



and an operator of arbitrary order. 

The problem appeared in connection with invariant construction 
of Darboux transformations for C. 



Mathematics Subject Classification (2010). Primary 70H06; Secondary 
34A26. 

Keywords. Joint differential invariants, gauge transformations, bivariate 
linear partial differential operator of arbitrary order. 



1. Introduction 

The present paper is devoted to Darboux transformations [TsaOO] for the 
Laplace operator, 

C = D x D y + aD x + bD y + c , (1.1) 

where a, b, c G K, where K is some differential field (see Sec. [2|. Operator C 
is transformed into operator C\ with the same principal symbol (see Sec. ^ 
by means of operator M. if there is a linear partial differential operator M 
such that 

NC = CiM. (1.2) 
In this case we shall say that there is a Darboux transformation for pair 
(£,jM); we also say that C admits a Darboux transformation generated by 
M. . We define the order of a Darboux transformation as the order of the M. 
corresponding to it. 
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Given some operator 1Z £ K[D] and an invertible function g £ K, the 
corresponding gauge transformation is denned as 

U^U 9 , V 9 = g~ l Rg . 

The principal symbol of an operator in i^[-D] is invariant under the gauge 
transformations. One can prove [Shel2] that if a Darboux transformation ex- 
ists for a pair (£,A4), then a Darboux transformation exists for the pair 
(£ 9 ,A4 9 ). Therefore, Darboux transformations can be considered for the 
equivalence classes of the pairs (C,M.). 

A function of the coefficients of an operator in K[D] and of the deriva- 
tives of the coefficients of the operator is called a differential invariant with 
respect to the gauge transformations if it is unaltered under the action of 
the gauge transformations. For several operators, we can consider joint dif- 
ferential invariants, which are functions of all their coefficients and of the 
derivatives of these coefficients. Differential invariants form a differential al- 
gebra over K. This algebra may be D-generated over K by some number 
of basis invariants. We say that these basis invariants form a generating set 
of invariants. For operators of the form known as Laplace invariants 

functions k — b y + ab — c and h = a x + ab — c form a generating set of 
invariants. 

In |SM08| we developed the regularized moving frames of Fels and 
Olver FQ99. f01vll[ [ManlOj for the individual linear partial differential op- 
erators of orders 2 and 3 on the plane under gauge transformations and 
obtained generating sets of invariants for those operators. In the present pa- 
per we extend those ideas and show that there is a finite generating set of 
invariants for the pairs (£,A4), where C is of the form (jl.ip and M. E K[D] 
is of arbitrary form and of order d. For M. of arbitrary order d but given in its 
normalized form without mixed derivatives, we find explicit general formulas 
for the basis invariants for the pairs (£,M) (Theorem 13. ip . 

The existence of such normalized forms for M. is implied by one of the 
theorems proved in |Shel2j . which can be re-formulated as follows. 

Theorem 1.1. Shel2] Let there be a Darboux transformation for pair (£, M), 
where C be of the form (1.1)) and M. £ of arbitrary form and order. 

Then there is a Darboux transformation for pair (£, M 1 ), where M! contains 
no mixed derivatives. 

2. Preliminaries 

Let K be a differential field of characteristic zero, equipped with commuting 
derivations d x ,d y . Let K[D] = K[D X , D y ] be the corresponding ring of linear 
partial differential operators over K, where D x , D y correspond to derivations 
d x ,d y . One can either assume field K to be differentially closed, in other 
words containing all the solutions of, in general nonlinear, Partial Differential 
Equations (PDEs) with coefficients in K, or simply assume that K contains 
the solutions of those PDEs that we encounter on the way. 
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Let / 6 K, and £ 6 if [I?]; by Cf we denote the composition of operator 
C with the operator of multiplication by a function /, while £(/) mean the 
application of operator C to /. The second lower index attached to a symbol 
denoting a function means the derivative of that function with respect to the 
variables listed there. For example, fi, X yy — d x d x d y f\. 

In the present paper we use Bell polynomials, 



B n ,k(xi,x 2 , . . . ,x n -k+i) 



E 



Jib'al • --in-fe+i! \ 1! 



Xi\0i ( x 2 \ 32 
2L 



Xn-k+l 

(n-k+l)l 



Jn-fc+1 



where the sum is over all sequences ji, J2, J3, ■ • ■ ,in-fe+l of non-negative in- 
tegers such that ji + J2 + • • ■ = k and ji + 2j 2 + 3j 3 + ■ ■ ■ = n. The sum 



B n (xx 



, • ■ • , 



) = B n<k (xi,X2, . . . ,X n -k + x) 



fc=l 



is called the n-th complete Bell polynomial, and also it has the following 
determinant representation: 



B n (xi,. . . , x n ) = dct 



(n—l\ (n—l\ (n — V\„ 

Xl { 1 )X 2 ( 2 )x 3 ( 3 )X A 

1 In— 2\ In— 2\ 

-1 Xi { 1 )X 2 { 2 JX 3 

n—3\ 







-1 



Xl 



( 1 )*2 



Xn 
Xn—1 

■ ■ x„_ 2 
T Xi 



3. Invariants for Normalized Darboux Transformations of 
Arbitrary Order 

Theorem 3.1. All joint differential invariant^ for the pairs (C,A4), where 

d 

C = D x D y + aD x + bD y + c and Ai = miD l x + m—iD y + too and where 

i=l 

mi £ K,i — —d, . . . , d and a, b, c £ K , can be generated by the following 
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2d + 3 basis invariants. 

m = a x -b y , 
h = ab — c + a x , 

Rj = J2^ ("W-rfC-ft, -d x (b), -%(b), , 

w=j ^ ' 

R -i = E m — ( W ) B w -j{~a, -d y (a), -d 2 y (a), . . . , -d^-\a)) , 

10=1 

+ m- w B w (-a, -d y (a), -dy{a), -dy~ 1 (a)) + m . 

Proof. We adopt the method of regularized moving frames |FQ991 IQlvlll 
IManlOj . Possible difficulties with the infinite dimensional case are addressed 
in OP08 . In this short paper we refer the reader to these works for the 
rigorous notation and for a justification of the method. 

For transformations £ i— > £ cx p( q ) j which implies the following group 
action on the coefficients of the operator, 

(a, b, c, ) -4 (a + aty, b + a x , c + aa x + ba y + a xy + a x a y ) , 

consider the prolonged action. 
Let us construct a frame 

P ■ (aj,bj,cj) (->■ g , 

at some regular point (x°,y°). Here a,j denotes the jet coefficients of a at 
(x° , y°), and they are to be regarded as the independent group parameters. A 
moving frame can be constructed through a normalization procedure based on 
a choice of a cross-section to the group orbits. Here we define a cross-section 
by normalization equations 

aj = , 
bx =0 , 

where here J is a string of the form x . . . xy . . . y, where y has to be present 
at least once, and there may be no x-s, and X is a string of the form x . . . x, 
where there can be no y-s. The normalization equations when solved for group 
parameters produces the moving frame section: 

aj = CLj-y , 

bx = b X -x , 

where J — y means that we take one y from the string J, and X — x means 
that we take one x from the string X. The first two fundamental differential 
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invariants can be then found 



by + a xy = by — a x = m = h — k 



p 

= c — ab — ab — a x + ab = c — a x — ab = h 



The remaining invariants of the generating set can be obtained using 
the constructed frame for the group acting on all the coefficients of the second 
operator in the pair, operator Ai since none of them has been used during 
the normalization process and construction of the frame. 

Consider a gauge transformation of each of the terms in the sum M. = 
Z)i=i miD x + m-iD l y + m : 

(mi£>J,) cxp ' a ' > = cxp(— a) ■ mi ■ D x o exp(a) , i ^ , (3-1) 
(m_ l Z?J / ) oxp(Q) = exp(-a) • • D l y o exp(a) , i ^ , 
(m ) cxp( - a) = m . 

Using the General Leibnitz rule |01v86] . for i ^ 0, we have 

(,n,I>;)"*'=exp(-o,. rai .t(-)^^>-^, 

then applying Faa di Bruno formula [S ROO] we continue 
,cxp(a) , , ^ (i\ \r^\ 9* exp(a) 



I — n V / +_i 



. . , da 1 

k=o v / t=i 



Bi-M {d x {a),dl{a), . . . , d x - k ~ t+1 (a)) ■ D k x 

i / - \ i—k 

E U ) E B i-UW, d 2 x (a),..., d x - k - t+1 (a)) ■ D k x 



k=0 



where Bi_k,t are Bell polynomials. Since only the terms Bi_k,t are summed 
with respect to t, we can rewrite the expression in terms of complete Bell 
polynomials: 

(m i B*)~ p(a) = m ([)B t - k (d x ( a ),dl(a), . . .,d x - k (a)) ■ D k . (3.2) 



fc=0 



Now we compute invariants Rj, j — 1, . . . , d from the statement of the the- 
orem as the coefficients at D J X , restricted on the constructed frame p. First, 
equality (|3.2j) implies that rtii appears only in Rj with j < i. Secondly, 
equality p. 21) implies that Rj must be a sum of the m, multiplied by some 
functional coefficients. The coefficient of rrii in Rj can be found from 
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by the substitution k = j. In this way, we can obtain the invariants 
R j =J2m w (™)B w ^(d x (a), d 2 x (a) : . . .,d^(a)) 



W=J 

d 



frame 



= E™« ( W )B w _ 3 (-b, -d x (b), -d 2 x {b), . . . , -B%-*-\b)) 

w=j ^ ^ ' 

for j = 1, . . . , d. Analogously, one can compute invariants R-j, j — 1, . . . , d 
from the statement of the theorem as the coefficients at D J y , restricted on the 
constructed frame p. 

We compute invariant Rq from the statement of the theorem as function 
A4(l) restricted on the constructed frame p. It has to be considered separately 
as this is the only invariant which contains both a and b. 



frame 



w—1 



+ Y. m AZ) Bw{d y (a) ' d v (a) ' ■ • ■ ' d v (a) } 



frame 



m 



= ^ m w B w {-b, -B x (b), -d 2 x {b), . . . , -d^-\b)) 

W — 1 

+ m- w B w (-a, -d v (a), -dy(a), . . . , -dy~ l (a)) + m . 



□ 



Theorem 3.2 (Alternative form of Theorem 13. All joint differential in- 

d 

variant^ for the pairs (A4,C), where A4 = ^""^ mjD x + m_j-D* + tuq and 

i=l 

C = D x D y + aD x + bD y + c, where rrn G K,i — —d, . . . , d and a,b,c G K , 
can be generated by the following 2d + 3 basis invariants. 



m 


= a x - 


&», 


h 


= ab — 


c + a 3 


Rq 


= m - 


d 

-£< 

i=l 


R j 


d-j 

=£( 

i=0 


\ 3 / 


R -j 


d-j 

=£( 

i=0 


+ i 
\ 3 / 



m i+j Pi(b) , j > 1, 
m_( i+ j)P;(a) , J > 1, 



2 with respect to gauge transformations 
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where 

Po(f) = 1 , 

p i (f) = -a i (f) ,ieN„, 

and the linear differential operator £1 is defined by 

UD.-m.it /-» 

U > \(D,-a)U) ,if } = a. 
Remark 3.3. Note the explicit forms taken by the first few operators Pf. 

W) = i, 

Pl(f) = -f, 

P 3 (/) = -/ M + 3//,-/ 3 , 

= -f xxx + i//,, - <i/,/ 2 + 3/ x 2 + .r , 

■:./•) - /, , , , + 5// xxx + 10/,,/, - 15/, 2 / - 10/,,/- + 10/,/'' - / 5 . 

5 

Example. Given the operator M = iriiD l x + m-iD y and an operator C in 

i=i 

the form with £ K,i = —5, . . . , 5, the following functions form a 

generating set of invariants. 

m = a x -b y , 

h = ab — c + a x , 

Ri = -2m 2 b + m x + (~3b x + 36 2 )m 3 + (-ib xx + Ylb x b - 46 3 )m 4 

+ (20b xx b + 156 2 ~ 306.6 2 - 5b xxx + 56 4 )m 5 , 

R 2 = m 2 - 3m 3 b + (-6b x + 66 2 )to 4 + (-10b xx - 106 3 + 306^6)777,5 , 

R 3 = 777,3 - 4m 4 6 + (-10b x + 106 2 )m 5 , 

i?4 = TO4 — 5777,56 , 
R5 = "7-5 , 

i?_ 5 = m_5 , 

i?_4 = 777_4 — 5m_5a , 

R-3 = 777 _3 — 4m_4a + (—10a,, + 10a 2 )m_5 , 

iJ_2 = 777,-2 — 3m_3a + (— 6a y + 6a 2 )77i_4 + (— lOa^y — 10a 3 + 30aj,a)777 _5 , 
= — 2777 _2<2 + 777_i + (— 3a y + 3a 2 ) 777 _3 + (— 4a yy + 12a y a — 4a 3 )m_4 
+ (20a TO a + 15a 2 — 30a y a 2 — 5a yyy + 5a 4 )m_5 . 

and finally 
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Rq =m - bmi - am-i + (-b x + b 2 )m 2 + (a 2 - a y )m- 2 + (3b x b - b a - b xx )m 3 
+ (3a y a - a yy - a 3 )m_ 3 + (6 4 + 3b 2 x - b xxx + Ab xx b - 66 x 6 2 )m 4 
+ (— 6a y a 2 + a 4 — a yyy + 3a 2 + &a yy a)m-4, 

+ {5b xxx b - b xxxx - b 5 - 10b xx b 2 + 10b x b 3 - 15b 2 x b + 10b xx b x )m 5 

+ (— 15a 2 a + lOdyyCLy + 5a y y y a — a y yyy — a 5 — 10a yy a 2 + 10aj,a 3 )m_5 . 
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